Robust quantization of a molecular motor motion in a stochastic environment. 
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We explore quantization of the response of a molecular motor to periodic modulation of 
control parameters. We formulate the Pumping-Quantization Theorem (PQT) that identi- 
^\ fies the conditions for robust integer quantized behavior of a periodically driven molecular 

machine. Implication of PQT on experiments with catenane molecules are discussed. 

(N 

Molecular motors are molecules capable of performing controlled mechanical motion. An ability 

to rotate its parts is a crucial function of a molecular motoi'^'^. It is challenging to control 

this motion in strongly fluctuating environment, experienced by any nanoscale system at room 

'~~| temperature. In the experimeniP, controlled rotation was implemented using 2- and 3-catenane 

1-1 

molecules that are made of interlocked polymer rings (n-catenane is made of n rings). Fig. [Ij 

^ shows geometry of a 3-catenane molecule and its 6 metastable states. Small (mobile) rings perform 

S transitions among 3 stations on the third, larger, ring. These transitions are caused by thermal 

2 fluctuations and, alone, may not lead to directed (clockwise or counterclockwise) motion on average, 
g 

I In experiments, the directed motion was induced by modulating the coupling strengths of mobile 

C rings to stations. This forced the smaller rings to orbit around the center of the larger ring while 

o 

remaining interlocked with it. 

In this Communication, we address an observatioiP'^'^ that a molecular motor can perform 
robust quantized operations, e.g. making a full mobile ring rotation per cycle in the control 
parameter space, even though the ring transitions are stochastic. We illustrate the phenomenon of 

m 
\6 

of flnding the number of rotations can be formulated in terms of stochastic motion on a graph 
X whose nodes and edges represent the metastable states and allowed transitions, respectively. 



integer quantization using a speciflc example of a 6-state stochastic model in Fig. [T} The problem 



• ^ The transition rates that satisfy the detailed balance can be written in Arrhenius form, i.e. they 

can be parameterized by the well depths Ei and potential barriers Wij = Wji, so that the kinetic 
transition rate from the node j to i is given by kij = e^^^^~^^^\ with /3 = l/ksT being the inverse 
temperature. Even if the rates can be written in Arrhenius form at any time, periodic changes of 
well depths and barriers result in a directed particle motion. The phenomenon is referred to as the 
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FIG. 1: Geometry and metastable states of a 3-catenane molecule. 

stochastic pump 

For the model in Fig. [Tj we further express Wij and Ei in terms of the directly controllable 
parameters, represented by the coupling energies e* of i-th mobile ring to j-th station, or the 
potential barriers between the station n, occupied by the z-th ring, and an empty station 
m. By requesting mobile rings be unable to occupy the same station simultaneously, and by 
comparing state energies and kinetic rates written in different parametrizations, e.g. Ei = e\ + e| 
or = e^i^i-^li) ^ we can relate the energies e*- and barriers to E^ and Wki in an 

effective 6-state model in Fig. [l| that is, E2 = e\ + el, E^ = £2 + el, E^ = + £\, E^ = £\ + el, 
Ee = el + el, Wu = el + w^g, W23 = + wi2> 1^34 = £2 + ^i3' ^45 = 4 + ^23^ ^^56 = 4 + ^i2> 
and Wie = e^ + ojI^. 

Consider a cyclic adiabatic evolution of e\ and eg- We introduce a pump current vector Q^, 
whose components Qf^ are average numbers of times the particle passed through links (i, j) during 
the cycle s of control parameters evolution, i.e. = J{t)dt, where J{t) is the instantaneous 
current vector and T is time of cyclic evolution of control parameters. The graph corresponding to 
Fig. |T] has one loop and hence Qf j+i = Q for any i, i.e. pump currents are the same for each link. 

Fig. [2] shows the dependence of the number Q of ring rotations per cycle on /?, obtained by 
solving the Master Equation in the adiabatic limit. It shows that, generally, the system's response 
to a periodic parameter variation is not quantized, however, in the low-temperature (3 ^ 00 limit, Q 
saturates to an integer value Q = 1. We have checked for a number of models that the phenomenon 
is generic. By choosing arbitrary closed contour in the space of control parameters, followed by 
choosing the remaining constant parameters randomly, integer response was always achieved in the 
(3 ^ 00 limit. 

The Pumping- Quantization Theorem rationalizes the observation and makes the following asser- 
tion: Consider any finite graph representing a Markov chain with kinetic rates written in Arrhenius 
form. If during the cyclic evolution of control parameters no degeneracy of the potential barriers 



3 




FIG. 2: Average number of rotations of mobile rings vs inverse temperature in a 3-catenane molecule 
after adiabatic evolution of control parameters along the contour e\ = —ef = cos(0), = —£2 = sin(0), 
e (0,27r); c^ig = 1/4, c^^^ = 1/2, c^^^ - -1/2, = 0, c^^g = -1/8, c^ig = -1/7, 4 = = 0. 

can encounter simultaneously with degeneracy of the minimal well depths, then in the adiabatic 
and after this low temperature (/? — > cxd) limits, the average number of particle transitions through 
any link of a graph per a driving cycle is an integer. 

We emphasize that by considering /3 — > 00 limit after the adiabatic approximation we assume 
that the particle has sufficient time to explore the whole phase space by making stochastic transi- 
tions before substantial change of control parameters can happen. Hence the number of rotations 
of the system per cycle is random and the quantization, that we discuss, appears only on average. 
Our arguments for the PQT are based on showing that various particle paths on a graph, that 
essentially contribute to the total current, are homologically equivalent to a single closed path, i.e. 
they differ from each other only by multiple transitions in both directions through some links. This 
is sufficient for a proof, since a single closed path can obviously pass through any link only an 
integer number of times. The proof of PQT includes three steps. 

(i) First, we prove the following identity 

(J, J^"""') = ef^^^^ Jij 4°°^ = 0, (1) 

{iJ} 

valid for the physical current J and any conserved current J"^""*^ that circulates in a loop of a graph 
and has equal values on any link of that loop. Summation in ([T]) runs over the graph links. Eq. ([T]) 
follows from a fact that the physical instantaneous current can always be represented in the form 
Jij = e~^^'^ {e^^^ 6pj — e^^'6pi), where 6pi is deviation of the probability on the i-th node from 
its equilibrium value. Taking any loop of a graph and summing {e^^^6pj — e^^^dpi) along its links 
gives zero, which is equivalent to ([T|. 

In the case of no barrier degeneracy, due to the e^^^^ factor in (fTl), the link {/c, /} with the largest 
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barrier inside any loop on a graph dominates the scalar product ([T| in the /? — > oo limit, and the 
only way to make this product zero is to conclude that lim^^oo('^A;/) = 0. Consider a time segment 
a without barrier degeneracies and let be the mean current integrated over the time of this 
segment. Integrating over time leads to inequality \Q%^ < X]{t l^^*'^*^"'^''''') \Q'ij\^ 

where summation is over all links of the loop except the link {A;,/} with the highest barrier, 
and max(...) means the maximum value of the expression during the given time segment. The 
factors max (e'^^^'J"^*^'^) become infinitely small in the /3 ^ oo limit but, for any given finite 
contour in the space of control parameters, IQ^ I remains finite even in the limit of infinitely long 
adiabatic evolution because pump currents, on average, depend only on a choice of a path of control 
parameters but do not depend explicitly on time of the evolutioii^. This leads to a stronger result 
lim/3_^oo IQfcil = 0. 

(ii) Here we note that the suppression of transitions through highest loop barriers during time 
interval without barrier degeneracies means that complex particle motion is restricted to a subgraph 

C X, referred to as the maximal spanning tree, that depends only on the ordering of non- 
degenerate barriers W. It is constructed from X step-by-step by eliminating the edge with the 
largest barrier that does not destroy the connectivity. Eventually, when none of the links can be 
removed without disconnecting the graph, we obtain the tree Xw. When energy Ej' approaches 
the lowest energy Ej and then becomes a new energy minimum, the particle travels from site j 
to site j' . Let ljji{Xw) be the unique shortest path that connects j to j' via Xw- Although 
transitions form j to j' are stochastic and can be done in a variety of ways, all paths from j to j' 
on Xw are homologically equivalent to ljj'{Xw)- 

(iii) To include events of barrier degeneracies, we partition the time period of a driving protocol 
into a set of small enough segments, so that for each segment we either encounter only the minimal 
well depth degeneracy or just the potential barrier degeneracy. We refer to them as to 0- and 1- 
segments, respectively. If necessary, we merge the consecutive segments of the same type to make 
the 0- and 1-segments alternating. In the /? ^ oo limit, no current is generated on the 1-segments, 
since the populations are concentrated in the node j with the lowest value of Ej but a nonzero on 
average pump current is possible only when the state probability vector changes^. The populations 
at the beginning and at the end of the 0-segment a are concentrated, respectively, in well-defined 
nodes ja and j^, determined by the neighboring 1-segments. Since the paths with nonvanishing 
probabilities from ja to j'^ belong to a tree X\y^, the total current Q", passing during the time of 
the 0-segment a has values 1 on any link on the shortest path lj^jr^{Xwa) and it has zero values 
on other links. Finally, the current per cycle = is generated by the concatenation of 
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FIG. 3: The 6-state model with triple degeneracy of barriers and double degeneracy of all well depths. 

the consecutive paths lj^j/^{X\Y^) that correspond to the a-segments. This exphcitly identifies the 
generated integer-valued current and completes the proof of the PQT. 

According to PQT, non-integer quantization in 3-catenane molecules is highly unlikely when all 
mobile rings and stations are different. In what follows we argue that robust fractional quantization 
can occur in systems, where, due to some symmetries, permanent (rigid) degeneracy of certain wells 
and/or barriers takes place. To illustrate the point we consider the 3-catenane system in Fig. [T| 
with a special symmetry: one of the mobile rings, called active, has three residence energies £j, 
available for control, whereas the other, passive, mobile ring has the same constant residing energy 
e = for all three stations. All barriers for both rings are constant and identical, equal to uj. The 
system is described by the Markov-chain model in Fig. |3] defined on the cyclic 6-node graph X with 
permanently degenerate wells E2j-i = E2j = £j, distinct barriers Wj = £j +uj with j = 1, 2, 3, and 
a permanently degenerate barrier Wq = to. The links {2j — l,2j} and {2j,2j + 1} describe the 
transitions when the passive and the active ring, respectively, switches to the non-occupied station. 

The total current can be calculated in a way similar to the integer-quantization case by 
partitioning the time period into a set of alternating 0- and 1- segments. Consider a segment a 
with no degeneracy among the parameters Wk with k = 0, 1, 2, 3. Similar to the integer-quantized 
case, the largest barriers still dominate the scalar product in Eq. ([T]), however, they can now be 
degenerate. Let Ap" be the change of the probability at site i during the 0-segment a and let 
maxa be the set of links {i, i + with the largest barriers during a. Then, in the low-temperature 
limit, Eq. ([T]) combined with the continuity equations leads to 

5^Q,Vi = 0' Ap'^ = Q'^_,^j - Ql^+„ j = l,...,5. (2) 

maxcj 

Eqs. (|2]) completely determine currents on each time segment. Note that solution of Eqs. ^ 
results in generally rational values for Q'jj^i because Ap" take values in a set (-1/2,0,1/2) and all 
other coefficients in equations are integers. 
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The protocols under consideration avoid the set Y of "bad" parameters, characterized by simul- 
taneous degeneracy of the lowest wells and highest barriers. For the model in Fig. |3j Y consists of 
4 lines Ya, a = 0, 1, 2, 3, in the space of e = (ei, £2, £3), as it is shown in Fig. [4| To understand the 
global phase diagram of the quantized response, we examine the currents Q^, the same for all links, 
that are generated by parameter motion along small contours enclosing Ya only one time. Finding 
the currents on all contributing 0-segments, as described above, and then summing them yields 
the currents of 1 and 1/3, as shown in Fig. |4] The independence of Qa on a variation of a contour 
around Ya allows the corresponding fluxes Qo = 1 and Qj = 1/3 for j = 1,2,3 to be associated 
with the lines. It also suggests a topological nature of the quantized current: the response to a 
general contour is given by the rational winding-index = Xla^aQaj where is the number 
of times the contour encloses the line Ya, taken with a proper sign depending on orientations (see 
Fig.|4]). 

There is an obvious analogy between topological properties of pump currents and the Aharonov- 
Bohm effect in quantum mechanics. In the latter, the phase of the electronic wavefunction changes 
upon enclosing a quasi-lD solenoid with a magnetic field by an amount proportional to the total 
flux of the field inside the solenoid. This analogy can be extended using the recent observation that 
stochastic pump effect is a geometric phenomenon^ ^ in a sense that integrated over time current 
can be written as a contour integral in the space of control parameters = <f^A- de, where A is a 
vector potential (gauge field) in the space of control parameters. This allows an effective "magnetic 
field" B = Ve X A to be introduced. Our explicit calculations show that, at low temperatures, 
the field B is localized in narrow tubes, carrying fluxes 1 and 1/3, which become Ya-hnes in the 
/3 ^ 00 limit. 

The arguments leading to rational quantization in our model with degeneracies can be applied 
to any graph with rigid degeneracy of some barriers and/or potential wells. For example, rigid 
degeneracies appear in the 3-catenane model with identical mobile rings. Similar considerations 
predict fractional quantization with a minimal ratio 1/2 for this system. We checked numeri- 
cally that this fractional quantization is robust when parameters are varied keeping mobile rings 
identical. However, it is destroyed as soon as mobile rings are made different. 

In conclusion, we have shown that current quantization in a stochastic system is a generic phe- 
nomenon and we identified the conditions for its observation. PQT directly applies to experiments 
with catenane molecules predicting integer-quantized response in a generic situation. We showed, 
however, that additional symmetries can lead to fractional quantization. For example, we predict 
1/2 quantization in a 3-catenane molecule that has two identical rings. We also showed that a 



FIG. 4: Topology of the degeneracy space and currents in the active/passive ring model 



3-catenane molecule can demonstrate fractional quantization with a minimal ratio of 1/3, which 
has not been observed previously. Quantization of a molecular motor response is topologically pro- 
tected. This robustness should have applications to the control of nanoscale systems, experiencing 
thermal fluctuations. 
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